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Here we show how the recent exact determination of the bond percolation threshold for the 
martini lattice can be used to provide approximations to the unsolved kagome and (3, 12 2 ) lattices. 
We present two different methods, one of which provides an approximation to the inhomogeneous 
kagome and (3, 12 2 ) bond problems, and the other gives estimates of p c for the homogeneous kagome 
(0.5244088...) and (3, 12 2 ) (0.7404212...) problems that respectively a gree with numerical results to 
five and six significant figures. 

PACS numbers: Ak 64.60 



Percolation Q, has provided some of the most in- 
triguing and difficult problems in statistical mechanics. 
Devised in 1957 by Broadbent and Hammersley it 
has served as the simplest example of a lattice process 
exhibiting a phase transition, and its study provides in- 
sight into more complicated physical models. 

The problem is very simply stated. Given any lattice, 
such as either of those shown in Fig. ^ we declare each 
bond to be in one of two states; open or closed. If a 
bond (although we could just as well consider sites) is 
open with probability p and closed with probability 1 — p, 
then clusters of various sizes will appear, with the average 
cluster size increasing as a function of p. In the limit 
of an infinite lattice there exists a critical value of this 
parameter, denoted p c and referred to as the percolation 
or critical threshold, where an infinite cluster will appear 
with probability 1. The value of p c is specific to each 
lattice. 

While the problem can be easily and precisely defined, 
exact solutions for thresholds (or anything else for that 
matter) have historically proved elusive, with results be- 
ing limited to a small set of lattices. Recent work 0, Q 
has significantly expanded this set, and in fact it was 
shown in Q that an infinite variety of problems are ex- 
actly solvable so long as their basic cells are contained 
between three vertices and are stacked in a particular 
self-dual way. Despite this recent progress, the most per- 
plexing unsolved problems still remain. In particular, the 
exact site percolation thresholds of the square and hon- 
eycomb (also called hexagonal) lattices, and the bond 
threshold of the kagome lattice are still unknown after 
nearly half a century of research in the field. The latter 
problem is one of the subjects of this Communication. 

The square, honeycomb, and kagome problems belong 
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FIG. 1: (a) the kagome lattice, (b) the (3, 12 2 ) (or 3-12) lat- 
tice. 



to an important subset of two dimensional lattices called 
the Archimedean lattices 6J , in which all sites are equiva- 
lent. There are 11 such graphs, and although both site Q 
and bond @ thresholds have been studied numerically for 
all of them, the only exactly solved problems are the bond 
thresholds of the square, honeycomb, and triangular [t| 
lattices, and the site thresholds of the triangular, kagome, 
and (3, 12 2 ) lattices. Note that finding the site threshold 
is a completely different problem from finding the bond 
threshold, and these last two site values are known only 
because of a trivial transformation from the honeycomb 
bond lattice — a transformation that does not help us in 
solving the bond problems. However, the (3, 12 2 ) lattice 
bears enough similarity to the kagome that the methods 
we present here will provide us with estimates for that 
bond threshold as well, one of which agrees with a recent 
numerical result 8] to its limit of precision, which is six 
significant figures. 

The bond threshold for the kagome lattice has pre- 
viously been the subject of several conjectures 0, 0, 
Il2l ITa ] . Using a method that predicted correct critical 
frontiers for the Potts model [Tj| on other lattices, Wu 
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[l5j conjectured that it would also work for the kagome, 
and, using the fact that percolation is the q — ► 1 limit 
of the Potts model [3, proposed that p c = 0.524430..., 
the solution of a polynomial we will encounter below. A 
few years afterward, and also in the context of the Potts 
model, Tsallis [Tl|, [jj} offered the competing conjecture 
p c = 0.522372..., employing an argument that also made 
correct predictions for other lattices. It was not until 
much later that both of these propositions were ruled 
out numerically 18] though fairly high precision was re- 
quired to exclude Wu's estimate. Tsallis also considered 
the (3, 12 2 ) lattice, and proposed p c (3, 12 2 ) = 0.739830... 

Aside from these various speculative methods, in which 
one makes conjectures that must be verified or rejected 
numerically, there are some rigorous results for the 
kagome and (3, 12 2 ) thresholds in the form of bounds 
on the values of p c . This work is largely carried out 
by Wierman and co-workers 19, 20], using a technique 
called substitution. The method is such that continual 
refinements are possible and the most current rigorous 
bounds are |2lj : 

0.522197 < p c (kagome) < 0.526873 , (1) 

and 

0.739773 <p c (3,12 2 ) < 0.741125 . (2) 

Various other quantities besides the standard perco- 
lation threshold have also been studied on the kagome 
lattice such as the mixed site-bond threshold [2^], a cor- 
related percolation threshold [2^, and an exact solution 
for the average cluster number on a kagome lattice strip 
[24) . among others. As already mentioned, the kagome 
Potts model has also received, and continues to receive, 
attention. In addition to the work already cited, some 
recent examples include |25|. and |26j] in which the con- 
jectures of Wu and Tsallis are discussed for various values 
of q. 

Here we show how a recent exact solution on a sim- 
ilar lattice, the martini lattice [Fig. Ufa)], can be used 
to provide precise estimates of the kagome and (3, 12 2 ) 
thresholds. 

The starting point of our analysis is the bond threshold 
for the martini lattice [Fig.[^a)]. For the general martini 
generator of Fig. |2[b), the method outlined in Q gives 
for the inhomogeneous critical surface 

1 - p\p 2 r 3 - p 2 p 3 r\ - P\P 3 r 2 - p\p 2 r\r 2 

- PiP 3 rir 3 - P2p3r 2 r3 + p\p 2 p 3 r 1 r 2 

+ P\P2P 3 r\r 3 + PiP2P 3 r 2 r 3 +PiP2rir 2 r 3 

+ PiP 3 rir 2 r 3 + P2Pd,r 1 r 2 r 3 - 2p 1 p 2 p 3 r 1 r 2 r 3 = (3) 

which was also reported recently in [Tof. Taking r» = 1, 
we get the result for the critical surface of the general 
honeycomb lattice |9j: 

1 - P1P2 - PiPi - P2PZ + P1P2P3 = , (4) 
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FIG. 2: a) The martini lattice, b) The assignment of proba- 
bilities for the inhomogeneous threshold 



and taking pi — 1 we get the following formula for the 
critical surface of the general triangular lattice |eJ : 

1 - n - r 2 - r 3 + rir 2 r 3 = . (5) 

For the first approach to the kagome lattice, we start 
with the inhomogeneous double-bond honeycomb lattice, 
whose unit cell is shown in Fig.^a). Replacing the bond 
with probability pi in the honeycomb lattice with a pair 
of bonds in series with probability PiU, we find from Q 
that the critical surface is given by 

l-piP2ht2~P2P 3 t2h-pip 3 t 1 t 3 +pip2p 3 t 1 t2h = . (6) 

Now consider the progression shown in Fig. El Starting 
with the double honeycomb lattice (a), changing every 
up star into a triangle gives the martini lattice (b), and 
changing the down stars gives the kagome lattice (c). 
The fact that the thresholds of the first two stages of 
this transformation are now known allows us to make 
guesses as to the way to reach the third. 

Comparing © with J3J) , it can be seen that the trans- 
formation 



tlt 2 - 


-> r 3 + rir 2 (l - r 3 ) 


{<) 


hh - 


-> r 2 + rir 3 (l - r 2 ) 


(8) 


t 2 h - 


-> n + r 2 r 3 (l - ri) 


(9) 


tlt2t 3 ~ 


r 1 r 2 r 3 + n.r 2 (l - r 3 ) 






+^2^3(1 - n) + rir 3 (l - r 2 ) 


(10) 



effectively turns the double honeycomb critical surface 
into the martini critical surface. These substitutions can 
be interpreted in terms of probabilities of connections 
between vertices on a triangle, i.e., tit 2 is the probability 
that a particular pair of vertices are connected on the 
star, and r 3 + r\r 2 (l — r 3 ) is the probability of the same 
thing on the triangle. The same transformations will 
also change the critical surface of the honeycomb lattice 
Q into that of the triangular (J5J — but note that we 
are not applying the star-triangle transformation here. 
In fact, these manipulations are largely formal, as the 
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FIG. 3: The transformation from the (a) double honeycomb, 
to the (b) martini, to the (c) kagome lattice. 



equation ltTU|l is not implied by (Q — ©■ Nevertheless, 
we conjecture that if we transform the down star the same 
way, we will be on the kagome critical surface. Using JJJ)- 
l|10|l with t i replaced by pi and by Si , we find that 
becomes 

1 - riS! - r 2 s 2 - r 3 s 3 - sxr 2 r 3 - s 2 rir 3 - s 3 rir 2 

- ris 2 s 3 - r 2 S!S 3 - r 3 sis 2 + sirir 2 r 3 + s 2 rir 2 r 3 
+ s 3 rir 2 r 3 + rir 2 SiS 3 + r 1 r 3 s 1 s 2 + r 2 r 3 sis 2 

+ r 2 r 3 sis 3 + nr 2 s 2 s 3 + r 2 si-s 2 s 3 + r 3 sis 2 s 3 

+ rir 3 s 2 s 3 + r 1 sis 2 s 3 - rir 2 r 3 sis 3 - rir 2 r 3 s 2 s 3 

- r 1 r 2 r 3 sis 2 - rir 2 sis 2 s 3 - rir 3 sis 2 s 3 

- r 2 r 3 sis 2 s 3 + rir 2 r 3 sis 2 s 3 = . (11) 

Setting all probabilities equal gives the condition 

1 - 3p 2 - 6p 3 + 12p 4 - 6p 5 + p e = , (12) 

with solution in [0,1] p c = 0.5244297175.... This result 
turns out to be identical to the conjecture made several 
years ago by Wu by different means. Subsequently, 
this value was found to be high numerically, but by only 
3 • 10" 5 Note that dJ is a plausible form for the 

kagome threshold: all the bonds are equivalent, setting 
any one probability to gives the correct threshold for 
the A lattice [the lattice that results when pi is set to 
1 in Fig. EJb)] , and setting all pi = 1 reduces the ex- 
pression to the triangular critical surface. It is difficult 
to imagine any other form that satisfies these conditions 
and remains linear in the probabilities, suggesting that 
the true general formula for the kagome lattice will not 
be linear in this way. 

The same procedure can also be used to find an approx- 
imate solution to the (3, 12 2 ) lattice. We start with the 
triple-bond honeycomb lattice, and transform the stars 
into triangles in the same manner as before (Fig. Q). 
There are nine probabilities in this case and the resulting 




FIG. 4: Progression from the triple-bond honeycomb to the 
(3, 12 2 ) lattice. 



inhomogeneous condition is 

1 - miTO 2 (r 3 + rir 2 - rir 2 r 3 )(s 3 + sis 2 - sis 2 s 3 ) 

- mim 3 (r 2 + rir 3 - rir 2 r 3 )(s 2 + sis 3 - SiS 2 s 3 ) 

- m 2 m 3 (ri + r 2 r 3 - rir 2 r 3 )(si + s 2 s 3 - Sis 2 s 3 ) 
+ mim 2 m 3 (rir 2 + rir 3 + r 2 r 3 - 2rir 2 r 3 ) 

x (sis 2 + S1S3 + s 2 s 3 - 2sxs 2 s 3 ) = . (13) 

Setting all m, = 1 gives (|llfl (in factored form), and 
setting all mi = m and ri = s, = r gives the equation 
for an inhomogeneous (3, 12 2 ) lattice with all triangle 
bonds having probability r and all linking bonds having 
probability m: 

1 - 3m 2 (r + r 2 - r 3 ) 2 + m 3 (3r 2 - 2r 3 ) 2 = . (14) 

Finally, letting r = m = p gives the equation for the 
homogeneous (3, 12 2 ) lattice, 

(l+p-2p 3 +p 4 )(l-p + p 2 +p 3 -7p 4 + 4p 5 ) = , (15) 

with solution on [0, 1] p c = 0.7404233179..., well within 
the bounds of @ . According to the numerical analysis of 
Parviainen 0, p c (3, 12 2 ) = 0.74042195(80). Our result 
is high by less than two standard deviations. Yet, we can 
get even better agreement with both of these results by 
taking a somewhat different route. 

In our second approach, we also compare the critical 
double honeycomb with the critical martini lattice, but 
we consider all bonds equivalent, in which case the dou- 
ble honeycomb threshold is po = \/l — 2sin7r/18 by ©■ 
Now, consider the martini lattice with p\ = p 2 = p 3 = p, 
and r\ = r 2 = r 3 = r. Equation implies that the 
critical surface is 

1 - 3p 2 (r + r 2 - r 3 ) + p 3 (3r 2 - 2r 3 ) = . (16) 

and taking p — po, we find that the critical value for r is 

r = 0.52440876529769 .. . . (17) 

That is, when one star with bond probabilities po is re- 
placed by a triangle with probabilities r, the system re- 
mains at a critical point (even though local correlations 
will necessarily be different because this is not a fixed 
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FIG. 5: The substitution of probabilities for the second 
(3, 12 2 ) lattice threshold estimate. 



point of the star-triangle transformation). If we conjec- 
ture that the system still remains at a critical point when 
we make the same replacement for the other triangle, 
then H17|l is an estimate for the p c of the kagome lat- 
tice. In fact, (|17|l is very close to the numerical result, 
p c = 0.5244053(3) 18], although outside the given error 
bars. 

It turns out that (|17|l is numerically identical to the 
value conjectured by Hori and Kitahara, which however 
is only available as a conference abstract without 
a derivation. Evidently, we have effectively duplicated 
the derivation of these authors. However, we can go fur- 
ther and use our argument to estimate the threshold for 
the (3, 12 2 ) lattice. Again we start off with the double 
honeycomb lattice at the uniform threshold of po, and 
compare to a critical martini lattice with p = po \fr [Fig. 
|Sfa)]. The argument works as in the kagome case, with 
the transformation to the (3, 12 2 ) lattice shown in Fig. 
Efb). The solution to yields 

r = 0.74042117858374... . (18) 



This result is within the error bars of Q and falls within 
the rigorous bounds of 0], which raises the possibility 
that the result is exact. Clearly, more precise numerical 
work for both lattices is called for. 

We can generalize our argument above for the in- 
homogeneous (3, 12 2 ) lattice with two probabilities m 
and r. The critical surface is determined by (|16|l with 
p = poy/m. When m = 1, this gives the kagome estimate 
(|17f) . when m = r it gives the homogeneous estimate 
(|18f) . and when r = 1 it gives the exact honeycomb re- 
sult m = Pq. The formula (|16|l (with p = Po^/m) can be 
compared with ijT4*|) , which though mathematically quite 
different, gives very similar numerical solutions. Finally, 
we note one last relation: if we require that the second 
terms of the two estimates (|14JI and (|16|l (which repre- 
sents two-point correlations) be the same, we get the sim- 
ple condition 

pl/m = r + r 2 - r 3 (19) 

which turns out to be identical to Tsallis' conjecture for 
this system. As mentioned above, however, the predic- 
tions of this formula are much farther from the numerical 
measurements than the predictions of i|14|) and i|16|l . 

In conclusion, we have shown that the results for the 
martini and honeycomb lattices can be used to make pre- 
cise estimates of bond percolation on the kagome and 
(3, 12 2 ) lattices, both long-standing problems in percola- 
tion theory. For the kagome lattice, we have reproduced 
the conjectures of both Wu and of Hori and Kitahara, 
while for the (3, 12 2 ) lattice we have apparently very pre- 
cise estimates. Perhaps these methods can point the way 
to finding rigorous thresholds for these lattices, and an- 
alyzing other unsolved lattices in percolation. 

The authors thank W. May and J. Wierman for pro- 
viding a copy of their preprint prior to publication. R.Z. 
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